Dynamical scaling for critical states: is Chalker's ansatz valid for strong fractality? 
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The dynamical scaling for statistics of critical multifractal eigenstates proposed by Chalker is 
analytically verified for the critical random matrix ensemble in the limit of strong multifractality 
controlled by the small parameter 6 1 . The power law behavior of the quantum return probability 
Pjv(r) as a function of the matrix size N or time r is confirmed in the limits t/N — > oo and 
N/t —¥ oo, respectively, and it is shown that the exponents characterizing these power laws are 
equal to each other up to the order b 2 . The corresponding analytical expression for the fractal 
dimension di is found. 
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Fifty years after the seminal paper by P.W. Anderson 
[l[ which predicted localization it is becoming a common 
wisdom that the multi-fractal scaling properties of criti- 
cal [H and near critical 0] states are not just an exotic 
theory exercise but an important physics reality 0, 
Since the pioneer's work by F.Wegner [fjj], it is known 
that at the Anderson transition point the local moments 
of the wave- functions VVi(r) scale with the system size L 
as (summation over n runs over the critical window near 
the mobility edge) 



in 



where d q is the fractal dimension corresponding to the 
q-th moment {multi- fractality). This implies Q that all 
scales of the amplitude I'i/v^r)! 2 ~ L~ a are present with 
the corresponding number of sites N a ~ L^ a \ The 
function f(a) known as spectrum of fractal dimensions 
is given by the Lcgendre transform of r(q) = d q (q — 1). 
One can distinguish between the weak multi- fractality 
where d q is close to the system dimensionality d and is 
almost linear in q for not too large moments q > 1, and 
the strong multi- fractality where d q <?C d for all q > 1. 

Recently the interest in Anderson localization has been 
shifted towards its effect for a many-body system of in- 
teracted particles. The simplest problems of this type are 
the " multi- fractal" superconductivity 0] and the Kondo 
effect For such problems the relevant quantity is the 
matrix element of (local) interaction which involves wave- 
functions of two different energies E m and E n : 

K(u, R) =Y, E^«( r + R )! 2 I'M 1 ")!" 8{E m -E n -u)). 



r n.rn 



(2) 

A promising playground to observe such effects at con- 
trollable strength of disorder and interaction are systems 



of cold atoms where one-dimensional Anderson localiza- 
tion has been already observed Q and the observation of 
their two- and three- dimensional counterparts is on the 
way. 

It was conjectured long ago by Chalker and Daniel 
[lfjl [Til and confirmed by numerous computer simula- 
tions [1, OH G3 that at E a > lo > A (A is the mean 
level separation): 

K(uj, 0) = C{u) cx {E q /ojY, fj, = l- d 2 /d. (3) 

The scaling relationship Eq. ([3]) can be viewed as a result 
of application of the dynamical scaling 



L — >• L u cx to d , 
to the earlier Wegner's result [l3[ 



K(0,R) cx (L/R) 



d—d 2 



(4) 
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with the simultaneous assumption that the indepen- 
dence is saturated for R < £; £ being the minimum scale 
of multi-fractality of the order of the elastic scattering 
length and Eq cx £~ d is the corresponding high-energy 
cutoff. This property is of great importance for electron 
interaction in the vicinity of the Anderson transition, as 
it leads to a dramatic enhancement of matrix elements 
compared to the case of absence of multi-fractality 0, 0] ■ 
Note that the dynamical scaling hypothesis Eq. Q im- 
plies that the energy scale u> corresponds to the length 
scale L u equal to the size of a sample where the mean 
level spacing is uj. However simple and natural is this 
hypothesis, it leads to a somewhat counter-intuitive con- 
sequence in the case of strong multi-fractality. Indeed, 
in the limit d% — > the exponent fj. in Eq.Q saturates 
at fi = 1. It signals of a strong overlap of two infinitely 
sparse fractal wavefunctions, while one would expect such 



2 



states not to overlap in a typical realization, similar to 
two localized wavcf unctions. At present there is no ana- 
lytical evidence of validity of Eq.Q in the limit of strong 
multi-fractality though numerical simulations seem to be 
very encouraging to it 0]. The main goal of this paper 
is to provide such an evidence. 

To achieve this objective we consider a model sys- 
tem described by a random matrix Hamiltonian with the 
multi- fractal eigenstates HH . The Gaussian ensemble 
of such random NxN matrices is defined by the variance 
of random (complex) entries 



= (\H nm \ ) - -j- 



2 U 



m 



2(l-e 



0. 



(6) 

and describes a long-range hoppings between sites of the 
one-dimensional lattice. The parameter b determines the 
strength of the fractality, with the weak multi-fractality 
at b S> 1 and the strong multi-fractality at b <C 1. It is 
the last limit that we will concentrate on in this paper. 
For convenience of the further calculations we introduced 
here a regularizing parameter e —¥ +0. 

Numerous computer simulations 0, 0, [l6| show that 
the above model gives an unexpectedly good quantitative 
description of statistics of the critical wave functions in 
the 3D Anderson model. Inspired by this fact we use the 
critical RMT, Eq.©, as the simplest and representative 
description of multi-fractality at the Anderson transition. 

To proceed further on, we note that the moment I2 
defined in Eq.([T]) and the dynamic correlation function 
C(lj) given by Eq.© arc related with the physical ob- 
servable known as the return probability: 



P N (t) 



dio 
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Indeed, plugging K(uj,0) defined by Eq.j2j) into Eq.© 
and taking into account that the average value of 



,-i(E n -E m )t 

level separation A <■ 
lim 

t/N^c 



in the limit t 



at a fixed mean 



- TV -1 we obtain at A > 1: 
P N (t) = h oc N- d */ d , 



(8) 



where in a particular case of the critical RMT the fractal 
dimension di = Y^n=i c m b m allows a regular expansion 
in the parameter b <C 1. The first term of this expansion 
ci was computed in Ref . [T^ . 

Alternatively, if A ~ A -1 tends to zero at a fixed large 
t, the power law in Eq.(j3|) is unbounded from below and 
one obtains at t ^> Eq 1 ~ b^ 1 : 



lim 



P N {t) oc t^ 1 oc t 



-d 2 /d 



(9) 



Thus the Chalker's ansatz Eq.© is equivalent to the 
identity: 



d In P N (t) 



dint 



N/t 
t —. 



d In P N (t) 



-> 00 

00 



9 In A 



t/N -> 00 
A -> 00 



(10) 



It implies that at large t and N the leading dependence 
of In Pjv(i) on t and N is logarithmic, and the coefficients 
in front of In t and In N are the same in both limits and 
equal to —0*2/ d. This is exactly what we are going to 
demonstrate in the present work. 

In order to reach this goal one needs to find the return 
probability at finite time t as well as its limiting value at 
t — > 00. In this Letter we calculate Pjv(i) using the virial 
expansion in the number of resonant states, each of them 
is localized at a certain site n. The virial expansion for- 
malism was developed in following the initial 
idea of Ref. 12 Oil . The supersymmetric version of the virial 
expansion [191 ] is formulated in terms of integrals over 
super-matrices. In particular, it allows us to represent 
Pjv(i) as an infinite series of integrals over an increasing 
number of super-matrices Q n associated with different 
sites n. As it was shown in Ref. fl9{ . the terms Pjy(t) in- 
volving integration over i different super-matrices result 
in the contribution to Pjv(t) = J2i ' CO °f the form: 

(t)=b i - 1 p%\bt), (11) 

with the known [l9[ explicit expressions for (bt) and 

It follows from Ea.([TT]) that 



lnP N (t)~bp%Xbt) + b 



p^(w)-i(p^(w)) ! 



(12) 



For Eq. (fTU|) to be valid, one requires that in the limit of 
large bt, N: 

(i) P P(bt) = - Cl ln(min{6i,A}) 

(bt)) 



(h) 



pf(bt) 



-C2 ln(min {bt, N}) 



(iii) terms proportional to In 2 (bt) and In 2 A cancel out 
in the combination (ii). 

This cancelation, as well as the logarithmic asymptotic 
behavior with equal coefficients in front of \n(bt) and In N 
in the two different limits, is not trivial. We show below 
that these properties arc hidden in the general structure 
of the virial expansion and in the power law dependence 
of the variance of the critical RMT v n ^o tx 1/n 2 , Eq©. 

We start by analyzing the explicit expression for 
P ( N ] (t) obtained in Ref.[l|: 



N , 

pP(t) = -2^J2\ v 

n=l 



te~ 




where erf(z) = /* e * d£. 

In the large t and N limit the sum over n in Eq. (flU)) is 
dominated by large n at any e > 0. Replacing the sum 
by an integral J2 n fo° dn one can represent p^ (bt) 
as a double integral: 



/o y jn- 1 n 



dy 



dn 
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where r = ^L, J" 2 (z) = \/2^z 2 (l - z 2 )e~ z , and e > 
ensures convergence at small y. Now we take the logarith- 
mic derivatives of p N (bt) as in Eq. pU|) and implement 
the limits N/t — > oo or t/N — > oo. The results are 

-2V2^t^ / ^(z 1 "'), -2^iV e / %T 2 (z), 
Jo z Jo z 

(15) 

respectively. Finally we take the limit e — > 0. One can 
see that both expressions in Eq. ([15[) coincide and 

ci=2V§W _J- 2 ( Z ) = _, (16) 
Jo z V 2 

provided that the operations of taking the limit e — > 
and integrating over z commute. 

Thus the validity of the Chalker's ansatz, in the form 
given by Eq. (|10[) , in the first order in b <C 1 is based on 
the symmetry of the integrand in Eq. (JTJJ) with respect to 



n and y at e = 0. However the necessary condition for 
this argument to work is the commutativity of integrating 
up to infinity and taking the limit e — > 0. Such a com- 
mutativity can be checked straightforwardly for Eq. ([15"]) . 

Now we proceed with b 2 contributions. Here the higher 
powers of In N or In r may arise and their cancelation in 
In Pjv(t) is of principal importance. The cancelation, if 
it happens, excludes in Eq.® presence of further addi- 
tive terms such as N~ kd2 or r~ kd2 with k > 1. Note 
that such terms, albeit not changing the leading large- 
TV, r behavior, could significantly alter the perturbativc 
expansion in b In N and 61nr. Their absence is a strong 
argument in favor of a pure power-law behavior Eqs. (|8I9|) . 

Using the representation for pff (bt) derived in Ref. flil ] 
one can cast the combination V^\bt) = p^ - (t) — 

5 (pa^( t )) hi the following form: 



V T =T~ 



V N = N 2t 



1 - e 
2e 

00 dp 



dxdyF 3 (x,y;e) - 



1 /l - e 



2 V e 



1 

/ ^i-2 £ k u^dxdyFaix^^e) [l - \x - y\fi K ) - - 




(l-\x\0«) 



(17) 



(18) 



where K = 1/(1 - e), V N = tim T /N^oc 7/v ( r ) and ^ = um w/r^oo P]v ( T )> and 

^ 3 (x,y;e) = J^da ^e a [h{Z x )h{Z y ) (f 3 - h){Z x - y ) + f 2 (Z x )f 3 (Z y ) (Af, - h)(Z x - y )] , (19) 

A(Z) = Z/VTTZ, h(Z) = 1/VT+Z, f 3 (Z) = Z/(l + Z)§, Z x = l/(a\x\ 2 ^). (20) 

I 



The contour C is the Hankel contour encompassing the 
negative part of the real axis $t(a) < 0. 

We present the cumbersome Eqs. (fl~7|) - (fT9| not only to 
give a flavor of real complexity of the calculations but also 
to uncover the ultimate reason for the Chalker's ansatz 
to hold. As in Eq. (jl5p . there is a certain similarity in 
Eqs. pTf and (fT5|) which can be traced back to (n — m)~ 2 
behavior of the variance Eq. © . To exploit this similarity, 
we use the approximate equality: 



'dp 



o P 



1-6 



/(o) 



df 

d/3 



Applying this formula to Eq. (jl8|) . one can see that the 
first term on the r.h.s. reproduces immediately Eq.(fT7|) 
up to the change of the pre-factor t 1 -* byN e . Like 
in the dimensional regularization calculus [2l| , after tak- 
ing the logarithmic derivatives w.r.t. r or N and taking 
the subsequent limit e — > 0, the terms oc e~ t 1 -" and 



cx e _1 N 2e tend to a constant determining the b 2 contri- 
bution to d 2 . 

The full calculation, however, is complicated by the 
presence of e~ 2 singularity in T 3 leading to the appear- 
ance of the additional contributions of the form e e~ 2 etc. 
Thus one has to keep not only the first term on the r.h.s. 
of Eq. (|2"T|) but the next two terms as well. An accurate 
account of all such terms shows that Eq.(fTU|) is valid to 
the b 2 order, and the fractal dimension d 2 is equal to: 



lib {nbf 



10 



56 
3\/3 



In 4 + 7T I 



0(b 3 ), (22) 



(21) 

where 



7T 



dipi dip 2 dip 3 



(COS ifi + cos f 2 ) (cos <f\ 
I can be evaluated numerically 



- COS If 2 + COS f 3 ) 



0.79426047250532455983. 



(23) 
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possibility of symmetric representation of the two differ- 
ent limits which can be traced back to the (m — n)~ 2 
dependence of the critical variance. 
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FIG. 1: Duality relation verified by numerical diagonalization 
of RMT Eq.©. The deviation from Eq.(J2U) is less than 1% 
for all values of B = (-Kb) 2* . 



However, we believe that this integral may have a geo- 
metrical meaning which is still evading our comprehen- 
sion and thus it can be evaluated analytically. It re- 
motely resembles the integrals arising in the problem of 
resistance of the regular 3-dimensional resistor lattice [22 1 
which arc known to have an intimate relation with the 
number theory 

Note that the leading terms in dzQ)) at small and large 
b Q can be represented in the form of the duality relation: 



d 2 {B)+d 2 (B- 



1, B = (tt&)24. 



(24) 



It appears that this relation is well fulfilled at all values 
of B (sec Fig.l). Our analytical result Eq. (|2"2l together 
with the result [23| which gives no 1/b 2 terms in d 2 (b) 
for large b, implies that the duality relation is not exact. 
However, its extremely accurate approximate validity is 
only possible because of the anomalously small value of 
the coefficient w 0.083 in front of (nb) 2 in Eq. (j22]) . 

In conclusion, we have shown that the Chalkcr's dy- 
namical scaling and its drastic consequence for strong 
correlations of the sparse multi-fractal wavefunctions is 
valid in the critical random matrix ensemble in the limit 
of strong multi-fractality di -C 1. We checked its valid- 
ity in the form of Eq. (|10[) up to the second order in the 
small parameter b that controls the strength of the multi- 
fractality. Specifically (i) we observed the cancelation of 
the In 2 N and In 2 r terms in In P/v(t) required by a pure 
power-law behavior; (ii) we demonstrated that the coef- 
ficients in front of the In N and In r terms are the same 
up to the order b 2 , and (iii) we found analytically the 
b 2 term in the fractal dimension d 2 . The validity of the 
Chanter's ansatz in the form Eq. pTj|) is encoded in the 
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